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Nuclear Quantum Effects 

Dongdong Kang,^ Huayang Sun/ Jiayu Dai/' * Zcngxiu Zhao/ Yong Hou/ Jiaolong Zcng/ and Jianmin Yuan^'^' ^ 

^Department of Physics, College of Science, National University of Defense Technology, 

Changsha 410073, Hunan, People's Republic of China 

'^ State Key Laboratory of High Performance Computing, 

National University of Defense Technology, Changsha 410073, Hunan, People's Republic of China 

(Dated: April 4, 2013) 

Nuclear quantum effects (NQEs) on the structures and transport properties of dense liquid hy- 
drogen at densities of 10~100 g/cm and temperatures of 0.1~1 eV are fully assessed using ab initio 
path-integral molecular dynamics simulations. With the inclusion of NQEs, ionic diffusions are 
strongly enhanced by the magnitude from 100% to 15% with increasing temperature, while electri- 
cal conductivities are significantly suppressed. The analyses of ionic structures and zero-point energy 
show also the importance of NQEs in these regime. The significant quantum delocalization of ions 
introduces expressively different scattering cross section between protons compared with classical 
particle treatments, which can explain the large alterability of transport behaviors. Furthermore, 
the energy, pressure, and isotope effects are also greatly infiuenced by NQEs. The complex behaviors 
show that NQEs can not be neglected for dense hydrogen even in the warm dense regime. 

PACS numbers: 62.50.-p, 66.10.-x, 72.15.Cz, 67.90.+Z 
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Dense liquid hydrogen plays a crucial role in under- 
standing the material behaviors under extreme condi- 
tions [1-3] and many applications in astrophysics and 
energy sources [4-8] . Ionic and electronic transport prop- 
erties of warm and hot dense hydrogen are key points in 
the dynamics of capsule implosion of incrtial confinement 
fusion (ICF) [4], modeling plasma processes [9] and the 
interior structure and evolution of giant planets and ex- 
oplanets [6, 8, 10], of which the core pressures are esti- 
mated even up to 19 Gbar [11], while the temperature 
may be down to thousands of Kelvin. Small change of 
transport properties will remarkably alter the dynam- 
ics of planet's evolution [12] and hydrodynamical pro- 
cesses in ICF [13]. Thanks to the recent progress in 
experimental techniques, one could get access to this 
previously inaccessible regime of phase diagram of hy- 
drogen through a laser-induced shock wave loading of 
precompressed samples [14, 15]. However, theoretically 
understanding the nature of hydrogen under the ultra- 
high pressures is still rare and a great challenge. Den- 
sity functional theory based molecular dynamics simula- 
tions, named quantum molecular dynamics or ab initio 
molecular dynamics (AIMD) is one of the most success- 
ful method for (warm) dense matter [16], but typically 
employ a classical-particle approximation for ions. Using 
AIMD, transport properties at high temperature such as 
ionic diffusion, electrical conductivity have been obtained 
with good accuracy [17-19]. Unfortunately, this approx- 
imation could lead to inaccurate results, especially when 
the nuclear quantum character is noticeable for light- 
atom systems [20-24]. 

The protons in condensed hydrogen easily exhibit 
quantum effects due to their low masses. These quan- 
tum effects include nuclear quantum fluctuations and 
exchange effects involving Bose-Einstein or Fermi-Dirac 
statistics [25] . Path-integral molecular dynamics (PIMD) 



[26-29] have provided a proper description of the nuclear 
quantum fluctuations beyond the harmonic level. Re- 
cently, a lot of simulations largely focus on the role of nu- 
clear quantum effects (NQEs) in the structures of dense 
hydrogen below several 1000 K, showing rich phases and 
dynamics [30-33] induced by quantum fluctuations or 
zero-point motion (ZPM). Compared with the classical- 
particle treatment, large deviation for ionic diffusion and 
thermal conductivity of para-hydrogen is found below 
32 K induced by NQEs [34]. People usually consider 
ions as distinguishable classical particles in warm dense 
regime (i.e., above 1000 K) since the quantum fluctu- 
ations or quantum collisions are generally regarded to 
be suppressed at high temperatures. However, even for 
the collision energies up to 10 eV, it has been shown that 
elastic collision cross sections for protons can be different 
significantly in the quantum theory [35] from the classical 
treatment. Nevertheless, NQEs are generally ignored in 
ab initio simulations of hydrogen in warm dense regime. 
This classical treatment for nuclear motion would be in- 
sufficient for dense hydrogen. In viewpoint of collision 
physics, both the ionic and electronic transport prop- 
erties depend strongly on the scattering cross sections 
between particles. In usual AIMD, the corresponding 
scattering cross sections are statistically obtained from 
classical particle dynamics simulations and therefore are 
expected to differ from the results with NQEs. 

In this Letter, we aim to fully assess NQEs on the 
structures and transport properties of warm dense hydro- 
gen using ab initio PIMD (AI-PIMD). The results show 
that the NQEs play an impressive role on both static and 
dynamic properties and can not be neglected for dense 
hydrogen even at the high temperature of 1 eV. 

Electronic structure calculations were carried out 
based on pscudopotential with plane-wave expansion 
of electronic wave functions [36, 37] and generalized- 
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FIG. 1. (color online). Temperature dependence of the self- 
diffusion coefficients of hydrogen atoms at 10 g/cm^. For 
comparisons, both the resuhs from AI-CMD and AIMD sim- 
ulations are presented. The inset is the Arrhenius plot of the 
self-diffusion coefficient. 



gradient approximation for exchange and correlation [38] . 
Electronic distribution is subjected to Fermi-Dirac statis- 
tics [39]. A periodic supercell including 250-432 atoms 
was employed according to different densities, which can 
ensure convergence of both ionic and electronic properties 
with good accuracy [4, 40]. A bcc lattice was used as the 
initial structure according to the results in Ref. 41. AI- 
PIMD calculations as well as their classical counterparts 
in this work were performed with modified Quantum- 
ESPRESSO package [42]. Within the framework of AI- 
PIMD, the structural and thermodynamical properties 
were calculated using the primitive scheme [27, 28], while 
the real-time quantum dynamics of nuclei was obtained 
through the ab initio centroid path-integral molecular 
dynamics (AI-CMD) [43, 44]. The self-diffusion coeffi- 
cients of ions were calculated from the slope of the mean 
square displacement obtained from AI-CMD simulations. 
Langevin thermostat [45] was employed to overcome the 
nonergodic problem, which not only produces a canonical 
ensemble and compensates the calculated errors [37] , but 
also gives us an efficient unified description from cold con- 
densed matter to hot dense regime [37, 46]. The Trotter 
number was set to 16 after a convergence test [47]. The 
proper time step and sufficient total steps were employed 
[48]. The electrical conductivity and optical absorbtion 
coefficient were calculated via the Kubo-Greenwood for- 
mula [49, 50] with ABINIT package [51]. 

The density of hydrogen under consideration ranges 
from 10 g/cm"^ to 100 g/cm^, corresponding to 0.3 < r^ < 
0.65 {rs denotes the Wigner-Seitz sphere radius). For 
each state in our simulations, the temperature ranging 
from 0.1 eV to 1 eV is well above the corresponding quan- 
tum degeneracy temperature [25, 52], thus the exchange 
of particles is negligible. 

Firstly, we calculated self-diffusion coefficients of hy- 
drogen atoms at 10 g/cur^ with temperatures from 0.1 



eV to 1 eV, as shown in Fig. 1. There are remarkable dif- 
ferences between the results from AIMD and AI-PIMD 
simulations. First of all, the self-diffusion coefficient at 
0.1 eV from the AI-PIMD simulations is more than one 
order of magnitude larger than the AIMD value. This 
amazing increment is likely to arise from the phase tran- 
sition from solid to liquid with the inclusion of NQEs. It 
is confirmed directly by ionic trajectory analyses (See 
the Supplemental Material for simulated movies [47]), 
where in the AI-PIMD simulations hydrogen exhibits ob- 
vious liquid character, whereas in the AIMD simulations 
hydrogen atoms remain bcc structure with large vibra- 
tional amplitudes. The difference introduced from NQEs 
is also clearly demonstrated by radial distribution func- 
tions (RDFs) of hydrogen nuclei in Fig. 2(a). We note 
that the RDF of hydrogen nuclei with AIMD method ex- 
hibits long-range ordered solid-like character at the low 
temperature of 0.1 eV, which is not visible in RDF ob- 
tained from the AI-PIMD calculations. In the quantum- 
mechanical perspective, it can be understood from that 
the protons tunnel through the energy barriers and move 
away from their equilibrium positions in bcc structure at 
0.1 eV, as appeared in hydrogen-bonded systems at low 
temperatures [20]. Our results accord with the conclu- 
sions in Ref. 41 that the melting temperature of hydrogen 
including NQEs is lowered compared with the classical 
particle treatment. 

Interestingly, the self-diffusion coefficients obtained via 
AI-CMD are substantially larger than the classical par- 
ticle value over the whole temperature range of 0.1-1 eV, 
as shown in Fig. 1. The difference between them is from 
100% at 0.3 eV, 38% at 0.5 eV to 32% at 0.7 eV. There 
is still a large difference of 15% even at the high temper- 
ature of 1 eV. It is not open-and-shut because the RDFs 
from the AIMD and AI-PIMD calculations are very close 
at 0.3 cV and almost identical when the temperature is 
up to 1 eV (see Fig. 2(a)). In addition, it is well-known 
that in the regime satisfied Boltzmann statistics, the dif- 
fusion is thermally activated and the temperature de- 
pendence of the self-diffusion coefficient D{T) generally 
obey the Arrhenius law, D(T) oc exp{-Ea/kBT) [53], 
where the activation energy Ea is the energy required for 
the hydrogen to surmount the potential barrier separat- 
ing the neighboring configurations. Here, the Arrhenius 
plot (see the inset in Fig. 1) shows that the self-diffusion 
coefficients from the AIMD calculations exhibit the Ar- 
rhenius behavior above 0.3 eV. The diffusion at 0.1 eV 
is smaller than the value according to the Arrhenius re- 
lation, arising from the solidification in the AIMD simu- 
lations as mentioned above. On the contrary, the results 
from the AI-CMD calculations are obviously larger than 
the estimates of the Arrhenius relation below 0.5 eV, in- 
dicating that quantum diffusion mechanism is becoming 
crucial in the temperature range. We note that the large- 
angle scattering between ions is dominant at such high 
densities considered here. The pronounced quantum nu- 
clear wave scattering calculation leads to a smaller large- 




FIG. 2. (color online). Radial distribution functions of hy- 
drogen nuclei from AI-PIMD (solid lines) and AIMD (dashed 
lines) simulations with different temperatures at 10 g/cm 
(a), and with different densities at 1 eV (b). The isotopes 
effects of deuterium (D) (long-dashed lines) and tritium (T) 
(dot-dashed lines) are also presented. The inset is the first 
peak of radial distribution functions at 100 g/cm'^ and 1 eV. 



angle scattering cross section between ions than the clas- 
sical particle scattering calculation, and thus increases 
the mean free path of ions. Therefore, the ions diffuse 
more easily even beyond the quantum tunneling regime 
and the classical particle treatment of protons substan- 
tially underestimate the ionic diffusion. In fact, when the 
indistinguishable effects are notable, much more signifi- 
cant effects will be visible even up to 10 eV [35]. This 
microscopic interpretation allows one to understand pro- 
foundly the transport properties of dense matter even 
at high temperatures where NQEs are conventionally ig- 
nored in the AIMD simulations. 

The largely reduced and broadened first peak of RDF 
from the AI-PIMD calculations at 0.1 eV and 10 g/cm^ 
indicates the significant nuclear quantum dclocalization, 
which becomes weaker with increasing temperature and 
can not be observed at 1 cV finally. In contrast, RDFs 
show much pronounced nuclear quantum character with 
increasing density (see Fig. 2(b)). When the density is 
increased up to 100 g/cm'^, there is a distinct difference of 
RDFs between the ATPIMD and AIMD calculations due 
to NQEs even at the high temperature of 1 eV. Mean- 
while, the isotope substitutions provide a clear demon- 
stration of isotope effects on nuclear spacial distribution 
of condensed hydrogen. Since the tritium nuclei is the 
heaviest one and has the shortest thermal de Broglie 
wavelength [47] , RDFs of the tritons are more structured 
than those of the other two isotopes, and more close to 
the results of the AIMD simulations for protons as the 
experimental observations [54]. 

Of our particular interest is whether the NQEs has 
evident effects on the electronic transport properties as 
profound as nuclei exhibits. To shed light on this issue. 
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FIG. 3. (color online). Electrical conductivity and opti- 
cal absorption coefficient from the AI-PIMD (solie lines) and 
AIMD (dashed lines) calculations for there state points, i.e., 
10 g/cm'' and 0.1 eV (upper panel), 10 g/cm^ and 0.3 eV 
(second panel), 10 g/cm^ and 1 eV (third panel), 100 g/cm'^ 
and 1 eV (lower panel). The error bars denote the standard 
deviation of averaging different atomic configurations. 



we calculated the electrical conductivity and optical ab- 
sorption coefficient of dense hydrogen. In order to avoid 
the artificial drop at low frequency due to the finite num- 
ber of atoms in the simulations, we employed the Drude 
formula to estimate the dc conductivity at zero frequency 
limit [16]. It can be clearly seen from Fig. 3 that with the 
inclusion of NQEs, the low-frequency dependent electri- 
cal conductivity (7[uj) exhibits different trend compared 
with the results from the AIMD simulations, which di- 
rectly leads to the different dc conductivity (7dc{^ -^ 
except for the state point of 10 g/cm^ and 1 eV. The dc 
conductivity from the AI-PIMD calculations arc largely 
suppressed compared with the AIMD value by 30% at 10 
g/cm^, 0.1 eV, 25% at 10 g/cm^, 0.3 eV and 18% at 100 
g/cm'^, 1 eV. Interestingly, a visible enhancement of op- 
tical absorption for hydrogen of 10 g/cm'^ at 0.3 eV and 
100 g/cm"^ at 1 eV, displayed in Fig. 3. It is well known 
that NQEs will introduce more disorder for the systems, 
lowering the melting point and critical point [33]. In fact, 
NQEs introduce ionic dclocalization and closer ionic dis- 
tances, as shown in Fig. 2, which produce more local- 
ized electrons around the nuclei [47, 55]. The stronger 
localized electrons might enhance the optical oscillator 
strength and optical absorption cross sections of rela- 
tive high frequencies. In addition, NQEs cause nonuni- 
form broaden of energy bands because of the anharmonic 
effects, resulting significant difference in the electronic 
density of states (DOS). In particular, the DOSs around 
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FIG. 4. (color online) . Differences of total energy between the 
AI-PIMD and AIMD calculations at a density of 10 g/cm^ (a) 
and a temperature of 1 eV (b). 



Fermi level decrease by NQEs [47]. 

To clarify NQEs on equation of state, we calculated 
the energies and pressures, and note that both the ener- 
gies and pressures including NQEs are higher than the 
values of AIMD simulations over the whole temperature 
and density range in this work [47]. While the ZPM and 
anharmonic effects account for the increase of ionic ki- 
netic energies and ionic pressures, the quantum spacial 
delocalization allows the nuclei to move closer and the 
Pauli exclusion principles enhances the kinetic energies 
of electrons that surround nuclei, thus the energies and 
pressures contributed by electrons rise inevitably. Fig- 
ure 4 shows the differences of total energies between the 
AI-PIMD and AIMD calculations. The energy difference 
^E is dropped at an exponential rate with increasing 
temperature, while raised linearly with increasing den- 
sity. We note that the zero-point energy (ZPE) is im- 
portant for determining the stable high-pressure struc- 
ture [56-58] and melting point of dense hydrogen [41], 
but precisely estimating ZPE beyond harmonic approx- 
imation is still a challenge in dense hydrogen [56]. The 
actual ZPE would be much larger than that estimated by 
harmonic treatment at high pressures due to the strong 
anharmonicity of nuclear motions [59]. Here we fitted 
the temperature-dependent energy differences shown in 
Fig. 4(a) to an exponential expression and obtained the 
ZPE at 10 g/cm^ is 2.24 eV/proton [60], which is much 
larger than that we could obtain from extrapolating the 
ZPE curve in Ref. 57 (about 1.78 eV/proton). It is 
due to that our PIMD calculations include anharmonic- 
ity of the nuclear motions naturally. In fact, ZPM can 
also strongly affect the optical absorption at low tem- 
peratures, inducing large corrections to the optical and 
electronic properties [61], which can explain partly the 
above differences of the electron and photon transport 
behaviors induced by NQEs. 

In conclusion, we have performed AIMD and AI-PIMD 
simulations to study the structure and transport proper- 
ties of dense liquid hydrogen. The complex transport be- 
haviors show that even when the NQEs have little effect 
on the structures (lOg/cm'^ and 1 eV), the ionic diffusions 
are also significantly enhanced due to the lower collision 



cross sections with the inclusion of NQEs. Electronic 
transport properties are less sensitive to NQEs here than 
ionic diffusions, but also seriously affected at relatively 
low temperature. The energy and pressure from the AI- 
PIMD simulations are higher than the AIMD value in the 
whole range of conditions here. The ZPE at 10 g/cm'^ for 
hydrogen is estimated by fitting the differences of total 
energies between the AI-PIMD and AIMD calculations 
at a temperature ranging from 0.1 to 1 eV. We can see 
that the quantum nuclear character will induce complex 
behaviors for both the ionic and electronic transports 
of dense hydrogen and can not be neglected even up to 
warm dense regime. 
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